July 6, 2012 



Review of D Semi-leptonic Decays 



Chunlei Lull] 

Department of Physics 
Carnegie Mellon University, Pittsburgh, PA, 15213 



Semi-leptonic D decays continue to play an important role in the field 
of flavor physics. During this presentation, recent measurements from 
pseudo-scalar to pseudo-scalar modes, pseudo-scalar to vector modes, and 
rare modes will be discussed. These results are important for many pur- 
poses, such as validating the machinery of lattice QCD, extracting CKM 
matrix elements, and searching for new physics and new interactions. 
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1 Semileptonic D Decays in the Big Picture 



One important task in the field of flavor physics is to over-constrain the CKM matrix. 
By doing so not only will we understand better the standard model (SM) physics, 
but we also have the opportunity to discover new physics. To best extract CKM 
matrix elements, we need inputs both from experimental data and from theoretical 
calculations. Fig. [1] shows our best knowledge about the CKM matrix up to today pQ. 
Our understanding about sin2/3 is the best thanks to clean theory and large statistics 
from experiments. In contrast, our understanding of both mixing and \V u f,\ have 
limitations from theoretical predictions. 
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Figure 1: Constraints in the (p, rj) plane. The red hashed region of the global 
combination corresponds to 68% C.L.. 

For example, \V u b\ is measured from B — > -nev which has the differential decay 
rate as the following: 

In order to extract the value of \V u b\ from the branching fraction measurement, we 
need information on the hadronic form factor f + (q 2 ). The latest world average of 
\V ub \ is P: 

|\4 6 | x 10 3 = 3.92 ± 0.09 ± 0.45, (2) 

where the uncertainties are from experiment and theory respectively. As we can 
see, the error due to theory is much larger than the experimental one. Our best 
understanding about f+(q 2 ) is from lattice QCD calculations. To validate improved 
lattice QCD calculations for the form factor in semi-leptonic B decays, we can use 
semi-leptonic D decays as a test and calibration. The differential decay rate function 
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for the semi-leptonic D decays is 



dT 

dq 2 




Vcs {d) \V KM \U(Q 2 )\ 



(3) 



24tt 3 



where a massless lepton is assumed. Since the CKM matrix elements |14 s (d)| can be 
obtained precisely from unitarity, we have a reliable method to check the lattice QCD 
calculations. 

The above example demonstrates how charm physics, in particular semi-leptonic 
D decays, can contribute to the larger picture of flavor physics. And more specifically, 
there are several advantages in using these decays: 

• Key modes have large branching fractions and are experimentally easily acces- 
sible at threshold experiments (where the missing neutrino can be inferred from 
missing four-momentum). 

• The theory of semi-leptonic decays is relatively clean, since the weak and strong 
interaction portions can be factored apart. 

• The pseudo-scalar to pseudo-scalar decay modes provide a simple and clean way 
to measure the form factor. And pseudo-scalar to vector decay modes provide 
access to even more form factors, if desired. 

• Rare semi-leptonic decays can also provide good paths to find new physics or 
new interactions. 



The experimental results of D semi-leptonic decays can be summarized in two main 
categories: the exclusive decays and the inclusive decays. Several recent measure- 
ments from the exclusive semi-leptonic D decays are discussed here. 

2.1 Measurements of D — > Kev and D — >• 7rev 

As previously discussed, these channels provide constraints on lattice QCD given 
CKM matrix elements or vice versa, and they have been measured in many exper- 
iments such as FOCUS, Belle, BaBar and CLEO-c. The main goal currently is to 
improve the precision for the Cabibbo-suppressed decay D — > -nev. 

The BESIII experiment has taken ~2.9 fb _1 ?/>(3770) data during the 2010 and 
2011 data runs. Using one-third of the data, a partially-blind analysis has been done 
with the D° — > Kev and D° — > -nev decays. Preliminary results were first presented 
at this Charm 2012 conference. 
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Recent Results 
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2.1.1 Event Selection 



The BESIII experiment takes data at the BEPCII symmetric electron-positron col- 
lider, with the -?/>(3770) produced at threshold. Using the double tag technique, several 
hadronic D decays are fully reconstructed at first. In case of multiple candidates, the 
candidate with minimum AE is chosen; AE is defined as 

AE = E can d — Ebeam, (4) 

where E can d is the energy of the reconstructed tag mode, and E^ eam is the beam 
energy. The following four hadronic D decays are used: D° — > K~n + , D° — > K~7t + tt , 
D° — > K~7t + 7t°7t° , D° — > K^7i + 7i^7T + . The rriBc distributions of these hadronic tags 
with our nominal AE requirements are shown in Fig. [21 where the beam constrained 
energy ijibc is defined as: 

m BC = \]El eam - \p\ 2 cand , (5) 
where p can d is the momentum of the reconstructed tag mode. 



mBC(D°-)Ki) mBC (tf-tKe 1 ) 




Figure 2: m BC distributions of D° K-tt + , D° R-tt+tv , D° K-tt+ttV, 
D° — > K~7t + 7t~7t + , with a AE cut of ±3<r for modes without ir° and — 4cr to +3<r for 
modes with 7r°. 

After hadronic D tags are found, msc are further required to be between 1.858 
GeV and 1.874 GeV before attempting to reconstruct signal candidates. The corre- 
sponding tag yields and tag efficiencies are listed in Tab. [TJ To search for the signal 
candidates, the following requirements are applied: 
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I mode 




Tag Yields 


Fraction(%) 


Tag efficiency 


£>° 


-)■ A'-TT- 


- 


159929±413 


20.7 


62.08±0.07 


D° 


-> A-7T- 


-7T° 


323348±667 


41.8 


33.56±0.03 


D° 


-> A-7T- 




78467±480 


10.1 


14.93±0.04 


D° 


— )■ A"~7T 


~7T~7T + 


211910±550 


27.4 


36.80±0.04 



Table 1: Tag yields and tag efficiencies using ~ 923 pb 1 of ^(3770) data from BESIII. 



• Only two good charged tracks with opposite charges left in the event. 

• One track is identified as an electron candidate, and the other track is identified 
as a kaon (pion) candidate for the K + e~v (7r + e~z/) mode. 

• The electron has the same charge as the kaon track from the tag side D. 

• The candidate is vetoed if the most energetic unmatched shower has energy 
greater than 250 MeV, in order to suppress backgrounds with extra tt°s. 

The energy and momentum of the missing neutrino is inferred by using: 

^rniss ^beara -^hadron -^electron (6) 

Pmiss ^tag Pfiadron P electron- (7) 

The number of signal events are obtained by fitting the U = E miss — P m i SS dis- 
tributions. Fig. |3] shows the U distributions and fit projections for the decay of 
D — > K + e~V and D — > 7r + e~z7. 




Figure 3: U distributions of D — > K + e z/(left) and D — > ix + e z/(right). Blue, 
green, and red curves are the total fit, signal fit, and background fit, respectively. 
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2.1.2 Measurements of Branching Fractions and Partial Decay Rates 



Given the signal yields obtained from fitting U distributions and signal efficiencies 
obtained from signal Monte Carlo, the branching fractions are obtained through the 
following equation: 



B 



N ohs 

sig 



my 



V N obs ' a f 



tag ^tag,sigl ^tag 



T obs,a 



where is the total number of signal yields with all tag modes combined, N°°g' 
is the observed tag yields for tag mode a, ef ag is the tag efficiency for mode a, and 
is the combined signal and tag efficiency for mode a. Preliminary results of 



-tag, sig 



branching fractions are listed in Tab. [2j 



mo 


de 


Measurement ( % ) 


PDG value (%) 


CLEO-c value (%) 


Tf 


— >■ K + e~V 


3.542±0.030±0.067 


3.55±0.04 


3.50±0.03±0.04 


D° 


— > n + e~V 


0.288±0.008±0.005 


0.289±0.008 


0.288±0.008±0.003 



Table 2: Branching fraction measurements using 
BESIII. 



923 pb^ 1 of ^(3770) data from 



In order to measure form factor, partial decay rates are measured in different 
q 2 bins, where q 2 is the invariant mass squared of the electron-neutrino system. 
D° — > K + e~V candidates are divided into nine q 2 bins (in GeV 2 ): [0.0,0.2), [0.2,0.4), 



[0.4,0.6), [0.6,0.8), [0.8,1.0), [1.0,1.2), [1.2,1.4), [1.4,1.6), [1.6, oo) , while D ->■ 7i + e~V 
are divided into seven q 2 bins: [0.0,0.3), [0.3,0.6), [0.6,0.9), [0.9,1.2), [1.2,1.5), [1.5,2.0), 
[2.0,oo). Signal yields in each q 2 bin are obtained by fitting U distributions in that q 2 
range. Using an efficiency matrix vs. q 2 , obtained from Monte-Carlo simulation, to 
correct for smearing, and combining with tag yields and tag efficiencies from previous 
studies, the partial decay rates are obtained, as shown in Fig. [H 

With measured partial decay rates, we can also obtain the f+(q 2 ) at the center of 
each q 2 bin i using: 



Mti) = ^VaJgW' (9) 

where |14<j(s)| can be taken from the PDG value, and pf is the effective value averaged 
over q 2 , calculated from: 

* I P 3 \U(Q 2 )\ 2 dq 2 nm 
Pl " f\U(d 2 )\ 2 dq 2 ' 1 ) 

where the three-parameter series parameterization for f+{q 2 ) is used from fitting data 
later. The f+(q 2 ) distributions are shown in Fig. [5] with theoretical curves overlaid [2]. 
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Figure 4: Partial decay rates measurement using individual tag modes (points) and 
all tag modes combined (histogram) for decay of D° — > K + e~V (left) and D° n^e~T7 
(right). 
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Figure 5: f+(q 2 ) distributions for the decay of D — > K + e~V (left) and D — > n + e~T7 
(right). Points are measured from data, red curves are the theoretical predictions, 
and blue curves represents variations within one statistical standard sigma. 



2.1.3 Measurements of the Form Factors 

Three different parameterizations of the form factor f+(q 2 ) are considered. 

The first parameterization, known as the simple pole model, is dominated by a 
single pole [3]: 



/+(0) 



l-q 2 /m 2 pole ' 



'IV. 



where the values of m po i e is predicted to be Mo* . Note that, however, that m po i e is 
floated in the fit in addition to /+(0) to improve fit quality. 

The second parameterization is known as the modified pole model [3]: 



U(q 2 ) 



/+(0) 



3 - 

pole 



l-a4 



pole 



(12) 



where m po i e is generally fixed to the D*^ mass. Both /+(0) and a are free parameters. 
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The third parameterization is known as the series expansion [I]. In the series 
expansion model, the branch cut in the complex q 2 plane is mapped onto a unit circle 
\z\ < 1, where 



z{q ,t ) = \ , (13) 



ft+ -q 2 + y/U=t Q 

where t± = (mp ± mp) 2 , and to is an Y rea l number smaller than t + . This transfor- 
mation expands the form factor about q 2 = to a s: 

1 00 k 

/+(?2) = PtfW.t.) £ "*<*>> I 2 * 92 '*"'] ' (14) 

where a; are real coefficients, P(q 2 ) = z(q 2 ,m 2 D «) for kaon final states, P(q 2 ) = 1 
for pion final states, and (p(q 2 ,to) is any function that is analytic outside a cut in 
the q 2 plane that lies long the x-axis from t + to 00. The standard choices are to = 
t+(l - (1 - t_/t+) 1 / 2 ) and, for 0: 



,,2., /% 2 ,0)\ 5/2 ( z(q 2 ,t ) Y 1/2 ( ztf,t.) \- 3/i t + -g 2 

0(g ' to) = v~ u^J li^j (^f 1 (15) 

Notice that the term has a limit of l/(4t + ) as q 2 goes to zero. 

With the above parameterization, the data is fitted by minimizing this \ 2 function: 

n 

X 2 = E (AT, - AG.;) Ctf (AT,- - AG,) , (16) 

where Ar, is the measured partial decay rate in q 2 bin i, AGi is the predicted partial 
decay rate in q 2 bin i, and C is the covariance matrix of the measured partial decay 
rates. The fitted results are listed in Tab. [3J 



2.2 Measurement of D + — > K 7r + e + u e 

The D + — > K^Ti + e + v decay includes contributions from different Ktt resonances, 
including the dominant K (892)° piece, as well as a non-resonant amplitude. The 
measurement of the K (892)° (lowest vector meson) contribution can be used to study 
hadronic transition form factors. Using 347.5 fb _1 of data recorded at the T(4S) the 
BaBar collaboration has analyzed this decay channel [5]. In addition to measuring 
the form factors, a comprehensive analysis of many other items is also performed. 
These include the properties of the K (892)° meson, the 5-wave contribution to the 
decay, the variations of the Kit S'-wave phase versus the Ktt mass, and a search for 
radially excited P-wave and .D-wave Kir resonances. 

The formalism of the differential decay partial width is taken from Ref. [6], which 
necessarily includes five variables: 
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Simple Pole 

D° ->■ K+e~v 
D — > 7r^e v 


/+(0)|V^)| 
0.729±0.005±0.007 
0.142±0.003±0.001 


ITT-pole 

1.943±0.025±0.003 
1.876±0.023±0.004 




Modified Pole 
D° ->• K+e-z/ 
L> — >■ 7r^e z/ 


/+(0)|Kd(.)| 
0.725±0.006±0.007 
0.140±0.003±0.002 


a 

0.265±0.045±0.006 
0.315±0.071±0.012 




2 par. series 
D° ->■ K+e-z/ 
L> — >■ 7r + e z/ 


f+mv cd (s)\ 

0.726±0.006±0.007 
0.140±0.004±0.002 


T\ 

-2.034±0.196±0.022 
-2.117±0.163±0.027 




3 par. series 
D° ->■ K+e-z/ 
L> — >■ 7r^e z/ 


f + mVc d (s)\ 
0.729±0.008±0.007 
0.144±0.005±0.002 


ri 

-2.179±0.355±0.053 
-2.728±0.482±0.076 


r-2 

4.539±8.927±1.103 
4.194±3.122±0.448 



Table 3: Fitter parameters from form factor measurements for D — > K + e v and 
D — > 7r^e v. 



• m 2 , the squared mass of the kaon-pion system; 

• q 2 , the squared mass of the electron-neutrino system; 

• cos(6k), where 9k is the angle between the kaon three-momentum in the kaon 
rest frame and the line of flight of the kaon-pion system in the D rest frame; 

• cos(# e ), where 9 e is the angle between the electron three-momentum in the 
electron-neutrino rest frame and the line of flight of the electron-neutrino system 
in the D rest frame; 

• the angle between the normals to the planes defined in the D rest frame by 
the Kix pair and the ev pair. 

2.2.1 Event Selection 

To reconstruct D + decays, all the charged and neutral particles are boosted to the 
center-of-mass system and the event thrust axis is determined. A plane perpendicular 
to this axis defines two hemispheres, and the candidate hemisphere consists of a 
positron, a charged kaon, and a charged pion. A vertex fit is performed and events 
with probability larger than are kept. To estimate the neutrino momentum, the 
K~it + e + v e system is constrained to the D + mass with estimates of the D + direction 
and neutrino energy obtained from all tracks and unmatched showers measured in the 
event. The neutrino energy is evaluated by subtracting from the hemisphere energy 
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the energy of reconstructed particles contained in that hemisphere. To further reject 
BB background and continuum background (arising mainly from charm particles), 
Fisher discriminants are created from relevant variables. This discriminant retains 
40% of signal events and rejects 94% of the remaining background. About 244 xlO 3 
signal events are selected with a final signal-to-background ratio of 2.3 . 

2.2.2 Measurement of Form Factors and K* (892)° Properties 

Form factors Fi,2,3 can be expanded into partial waves to show explicit dependence 
on 6k- Considering only S, P and D waves, this gives: 

F l = Fio + F u cos9 K + F 12 3C ° S 9 * - 1 (17) 

1 [3 
F 2 = ~7|^2i + y ^F 22 cos6k 

1 [3 
F 3 = —^=^31 + y -F 32 cos6 K 

where Fi characterizes the 5-wave contribution, and Fn and F i2 correspond to the 
P— and D-waves, respectively. It is also possible to relate these form factors with 
the helicity form factors H ±, which can be in turn related to the two axial- vector 
form factors Ai^ 2 (q 2 ) and the vector form factor V(q 2 ). For the q 2 dependence, the 
simple pole parameterization is used: 



V{q 2 ) 

AM 2 ) 

A 2 {q 2 ) 



V{0) 



1 — q 2 /rriy ' 

1 — q 2 /m 2 A ' 
A 2 (Q) 

2 ' 



q 2 /m A 



where my and are expected to be close to mo. = 2.1 GeV and m^j = 2.5 
GeV, respectively. In the analysis, ratios of these form factors, evaluated at q 2 = 0, 
r v — V(0)/Ai(0) and r 2 = A 2 (0)/A 1 (0) are measured by studying the variation of 
partial decay rate versus kinematic variables; m A is allowed to fit while my is fixed. 
For the mass dependence, in case of the K*(892), a Breit-Wigner distribution is used. 
Tab. H] shows the fit results considering three different models: 

• S + iT(892)°, a signal made of the 7T(892)° and S -wave components. 

• S + iT(892)° + K*(1410)°, a signal made of the iT(892)°, 7T(1410)°, and 
S-w&ve components. 
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• S + K*(892)° + F*(1410)° + D, a signal made of the if*(892)°, K*(1410)°, S- 
and D-w&ve components. 

The fractions of signal components in different models are measured and shown in 
Tab. [5j The second model is considered as the nominal fit to data, with systematics 
error also listed in the table. 



Variable 


S + K* (892)° 


S + K* (892)° 
iT(1410)° 


S + K* (892)° 
K*{U10)° + D 


m K * {892 ) (MeV) 


894.77±0.08 


895.4±0.2±0.2 


895.27±0.21 


r^ (892) (MeV) 


45.78±0.23 


46.5±0.3±0.2 


46.38±0.26 


r BW (GeV 1 ) 


3.71±0.22 


2.1±0.5±0.5 


2.31±0.20 


m A (GeV) 


2.65±0.10 


2.63±0.10±0.13 


2.58±0.09 


rv 


1.458±0.016 


1.463±0.017±0.031 


1.471±0.016 


r 2 


0.804±0.020 


0.801±0.020±0.020 


0.786±0.020 



Table 4: Values of form factor and K (892)° parameters with different models. 



Component 


S + K*(892)° 


S + K* (892)° 
K*{U10)° 


S + K* (892)° 
K*{U10)° + D 


S'-wave 


5.62±0.14±0.13 


5.79±0.16±0.15 


5.69±0.16±0.15 


P-wave 


94.38 


94.21 


94.12 


iT(892)° 
7T(1410)° 


94.38 


94.11±0.74±0.75 


94.41±0.15±0.20 





0.33±0.13±0.19 


0.16±0.08±0.14 


D-wave 








0.19±0.09±0.09 



Table 5: Fractions (in percent) of signal components with different models. 



2.2.3 Measurement of S-wacve Phase 

Using the nominal signal model containing S'-wave, i^*(892)° and K (1410)° compo- 
nents, the phase in the mass- dependent S'-wave amplitude is measured and compared 
to LASS measurements. The S-wave phase is assumed to be a constant within each 
of the considered Kir mass intervals. The measurement of phase variation is shown 
in Fig. O 

2.3 Measurement of D°/D + pev and D + — > ujev 

Using 818 pb _1 of data taken at the ?/>(3770), the CLEO-c collaboration has measured 
branching fraction and form factor for the decays of D°/D + — > pev and the branching 
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Figure 6: Points (crosses) give the S'-wave phase variation using the nominal model. 
The phase is assumed to be constant within each mass bin, and parameters of the 
K (892)° are fixed to values from the nominal fit to data. The solid line corresponds 
to the parameterized S'-wave phase variation in the nominal fit. LASS data are shown 
in blue. 

fraction for D + — > uev [7J. The precision of the branching fractions is improved, and 
the form factor result of the D°/D + — > pev is the first measurement on the Cabibbo- 
suppressed pseudo-scalar meson to vector meson transition in semi-leptonic D decay. 

2.3.1 Event Selection 

The double-tag technique is used in this analysis by reconstructing a D tag in the fol- 
lowing hadronic final states: K + tt~ , K + tt^7T , and K + 7t~7t~tt + for neutral tags, and 
K®7T~ , K + 7t~tt~ , K^n^n , K + 7t~7t~7t , K®7r~ 7i~ 7i + , and K~K +r R~ for charged tags. 
In case of multiple candidates in the same tag mode, the candidate with minimum AE 
is chosen. Once a tag is identified, certain AE and vn^c cuts are required. The un- 
used tracks and showers are then searched for a candidate e + along with a p - (7r~7r°), 
p°(7r + 7r~), or u;(7r + 7r~7r ). The p candidate is required to have invariant mass within 
150 MeV from the nominal PDG mass. The combined tag and semi-leptonic candi- 
dates must account for all reconstructed tracks and unmatched showers in the event. 
To remove multiple candidates in each semi-leptonic mode, one combination is chosen 
per tag mode per tag charge, based on the proximity of the invariant masses of the 
p°, p + or uj candidates to their PDG masses. The number of semi-leptonic decays 
are obtained by fitting the U = E miss — \p m iss\ distributions, where E miss and p m iss 
are energy and momentum of the missing neutrino, and they can be inferred from all 
other measured particles. The U distributions are shown in Fig. [71 
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Figure 7: U distributions of (a) D° — > p~e + v ei (b) D + — > p e + v e and (c) D + — > ue + i / . 



2.3.2 Measurements of Branching Fraction 

To measure the absolute branching fraction, Monte-Carlo samples are used to deter- 
mine tag and signal efficiencies. Together with the tag and signal yields, the same 
equation as Eq. [8] can be used to give the branching fraction. The measured results 
are listed in Tab. |6j 



Decay mode 


e(%) 


1 ' sig 


Bsl 


J B 5L (ISGW2) 


B S l(FK) 


D° p-e + u e 


26.03±0.02 


304.6±20.9 


1.77±0.12±0.10 


1.0 


2.0 


D+ p°e+z/ e 


42.84±0.03 


447.4±24.5 


2.17±0.12i°i| 


1.3 


2.5 


D + ->• ue + v e 


14.67±0.03 


128.5±12.6 


1.82±0.18±0.07 


1.3 


2.5 



Table 6: Branching fractions for D° — > p~e + v e , D + — > p°e + z/ e and D + — > we + z/, 
from the CLEO-c analysis and two model predictions: ISGW2[S] and FK[5|. The 
uncertainties for signal efficiency e and signal yields N sig are statistical only. The 
efficiency includes the p and u branching fractions from the PDG. 



2.3.3 Measurements of Form Factor 

A form factor analysis is performed for D°/D + — > p~ / p°e + v e decays. The mechanism 
of this decay is similar to the D + — > K-nev discussed in Sec. I2.2[ except only P-wave 
is considered in this case. Three dominant form factors, two axial and one vector, A\, 
A 2 , and V, are used to describe the hadronic current. A simple pole model is assumed 
with the pole mass fixed as M^*^-) = 2.01 GeV and M^*(i+) = 2.42 GeV for the 
vector and axial form factors, respectively. A four- dimensional maximum likelihood 
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fit [10] is performed in the space of q 2 . cos9 n , cos8 e and Xi an d a simultaneous 
fit is made to the isospin-conjugate modes D° — > p°e + i> e and D + — > p~e + v e . The 
ratios of form factors evaluated at q 2 = 0, ry = ~x^jyj an d r% = ^jjjy are obtained: 
rv = 1.48 ± 0.15 ± 0.05 and r 2 = 0.83 ± 0.11 ± 0.04. Using \V cd \ = 0.2252 ± 0.0007, 
r D o = (410.1 ± 1.5) x 10- 15 s, and t d + = (1040 ± 7) x 10~ 15 s, from PDG 2010, 
form factor ratios and branching fraction results are combined to obtain : Ai(0) = 
0.56 ± 0.01±8;8|, A 2 (0) = 0.47 ± 0.06 ± 0.04, and 1/(0) = 0.84 ± 0.09±g;g|. The fit 
projection is shown in Fig. |8l where the difference between data and the fit projection 
for cosOjr may be due to S"-wave interference. 




cos e e % 

Figure 8: Projections of the combined p~ and p° data (points with statistics error 
bar) and the fit (solid histogram). The dashed lines show the background distribu- 
tions. 

2.4 Search for the Decay of — > uje + v 

Using data collected at a center-of-mass energy yfs = 4170 MeV, CLEO-c searched 
for the decay of Df — > ue + u, which can probe the four-quark content of the [TTj . 
The integrated luminosity of the data used in this analysis is 586 pb~ x , or 0.6xl0 6 
DfD~ inclusive decays. Recent work by Gronau and Rosner [12J concludes that any 
value of the branching fraction for — > ue + v exceeding 2xl0 -4 is unlikely to be 
from bj — cf) mixing. 

2.4.1 Event Selection 

About 95% of the e + e~ — > D^X events are formed from the following exclusive 
reactions [T3] : 

e+e-^ D+D*-, e+e~ -+ DjD*+, (19) 
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with equal amounts of them produced, and about 95% of the D* decays through D s ^y, 
To select signal events, the double tag technique is used. Eight D s tag modes 
are reconstructed: K®K~, K + K~tt~, K*~K*°, tt + 7i~7i~ , T]tt~, T]p~, 7r~?/(?77r + 7r~), 
7r~?/(p7). The four-momentum of a tag candidate is defined by (E tag ,p tag ), and the 
recoil mass M rec is defined as 

M rec = ^/{E b -E tag y-(p b -p tag y, (20) 

where (E b ,p b ) is the four-momentum of the colliding beams. The M rec distribution 
peaks at zero only for events where the photon is associated with signal side, but 
even if the photon is associated with the tag side, M rec is kinematically constrained 
such that \M rec — Md*\ < 55 MeV. So only events passing with this wide cut are 
retained. The distribution of invariant mass of the tag M tag is fitted with a double 
Gaussian function, and the weighted resolution a is used to define signal regions. All 
tag modes are required to to be within 2.5a from the peak position except for the 
rjp mode where it is selected with 2a. A second kinematic constraint is used, on the 
MM* 2 defined as 

MM* 2 = (E b — E tag — E^) 2 — (p b — p tag — p 7 ) 2 , (21) 

where MM* 2 is expected to peak at M^ s (independent of whether the 7 is from the tag 
or signal side). The selection criteria on MM* 2 is determined from a two-dimensional 
binned likelihood fit the (MM* 2 ,M tag ) space. 

The signal is then selected by requiring one positron candidate with charge op- 
posite to the tag charge, two charged pions of opposite charges, and no extra good 
tracks. A yet-unused good 7r° candidate is also required. There can be multiple 7r° 
candidates; if so, the one with lowest \ 2 — [(-^77 — M^/a^] 2 is chosen, where 
M 77 is the photon-photon mass, and a 77 is the calculated resolution. The positron, 
charged pions, and n° are added together to get a four- vector (E s ,p s ). With informa- 
tion from beams and tag side, the missing-mass squared of the neutrino is obtained: 

MM 2 = (E b -E tag -E 7 -E s ) 2 -(p b -p tag -p 7 -p s ) 2 . (22) 
2.4.2 Final Fit and Results 

To extract the signal yields, the mass of the 7r + 7r~7r° system, M3, is used with a cut on 
MM 2 as -0.05 < MM 2 < 0.05 GeV 2 . The M 3 distribution is shown in Fig. El where 
one can see peaks for the r\ and 0, but no clear signal peak at the u. An upper limit 
on the branching fraction of B(Df — > ue + v) < 0.20% is obtained at 90% confidence 
level. 
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Figure 9: M 3 (M(3tt) distribution. Solid line: signal selection, after M tag sideband 
subtraction. Dotted line: M tag sideband contribution. 



3 Summary 

Experimental measurements of semi-leptonic D decays have been studied very suc- 
cessful during the past few years, with contributions from experiments such as FO- 
CUS, Belle, BaBar and CLEO-c. During this presentation, recent results on semi- 
leptonic D decays have been discussed, with topic covering several topical aspects, 
i.e., pseudo-scalar to pseudo-scalar modes, pseudo-scalar to vector modes, and rare 
modes. 

Since the start of running in 2008, the newest of these experiments, BESIII, has 
taken about 2.9 fb^ 1 of data at ?/>(3770). With peak luminosity reaching more than 
6 x 10 32 (60% of the designed luminosity), BESIII is poised to take more data at 
■^(3770) and in the higher D s energy region. Using part of the data, BESIII has 
presented preliminary results of the D° — > K/itev decays. Results from the full 
dataset and other modes are coming in the near future. 
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